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LONG-TERM  GOAL 

The  ultimate  goal  of  our  research  is  the  development  of  a  complete  ice  dynamics  model  that  will 
include  lead  direction  and  ice  thickness  distribution  in  refrozen  leads. 

OBJECTIVES 

Northwest  Research  Associates,  Inc.  (NWRA)  objectives  during  this  past  year  were  to  conduct  the 
necessary  data  reduction  and  analysis  of  stress  data  to  check  the  new  model  and  to  share  the  data  with 
Dr.  R.  S.  Pritchard  for  his  use  in  model  development. 

APPROACH 

Our  approach  is  to  use  measured  sea  ice  stress  and  motion  to  study  the  behavior  of  refrozen  leads. 

WORK  COMPLETED 

All  of  the  SIMI  stress  data  has  been  processed  to  give  geophysical  stress  resultants  for  comparison  with 
the  ice  model.  One  paper  has  been  prepared  and  presented  at  the  IAHR  14th  International  Symposium 
on  Ice,  in  Potsdam,  NY. 

RESULTS 

NWRA's  research  shows  that  in  the  interpretation  of  sensor  oil  pressure  as  ice  stress,  the  two  most 
significant  corrections  are  for  the  time-dependent  behavior  of  the  ice  and  the  thermally-induced  stress. 
Other  possible  data  errors,  such  as  electrical  drift  correction  for  gauge  temperature  variation  in  the 
gauge  stiffness,  were  all  investigated  and  shown  to  be  minor  in  comparison  to  the  effects  mentioned. 

Thin,  circular,  oil-filled  flat-jack  ice-stress  sensors  were  used  to  measure  sea  ice  stress  during  ONR's 
SIMI  Program  (Coon  et  al.,  1995a;  Coon  et  al.,  1995b).  These  sensors  were  calibrated  in-situ  in 
compression  tests  at  two  speeds;  rapid  loading  (around  100  kPa/min.),  which  produced  a  square-wave 
pressure  history;  and  slow  loading  (around  one  kPa/min.),  which  produced  a  ramp  pressure  history. 

The  load  rate  of  the  ramp  tests  more  closely  matched  geophysical  load  rates,  while  the  rapid  load  rate 
of  the  square-wave  test  is  fast  enough  to  measure  the  elastic  response  of  the  ice-sensor  system.  Based 
on  the  test  results,  it  was  decided  to  model  the  response  of  the  sensor  to  ice  stress  with  a  dynamic  or 
time-dependent  type  inclusion  factor  of  the  form 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 
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where  a  is  the  ice  stress  normal  to  the  sensor  or  the  air  pressure  in  the  air  jack,  P  is  the  oil  pressure  in 
the  sensor  (with  pressure  position),  and  A,  B,  and  x0  are  constants.  Elastic  inclusion  for  flat  jack 
gauges  for  ice  are  presented  by  Chen  (1981).  Hamza  and  Blanchet  (1984)  used  an  incremental  analysis 
for  the  viscoelastic  bridging  response  of  a  pressure  sensor  embedded  in  an  ice  sheet.  Their  approach  is 
similar  to  what  is  being  done  here. 

By  plotting  the  air  pressure  versus  the  oil  pressure  found  in  the  square-wave  tests,  the  effective  elastic- 
inclusion  factor,  A,  was  obtained  for  the  sensor  installed  first-year  sea  ice.  In  many  of  these  plots,  a 
change  in  slope  occurs  at  10  to  20  kPa  oil  pressure.  We  judged  that  the  stress  sensor  readings  are 
unreliable  whenever  the  oil  pressure  is  less  than  20  kPa  and,  therefore,  concluded  that  these  sensors  do 
not  measure  tensile  stresses  or  small  compressive  stresses.  The  creep  terms  in  the  dynamic  inclusion 
factor,  B  and  t0,  were  determined  using  the  linear  up-ramp  portion  of  each  ramp  test  along  with  the 
effective  elastic  inclusion  factor  from  a  prior  square-wave  test  at  the  same  installation.  The  measured 
air  pressure  for  the  whole  test  was  compared  with  the  air  pressure  calculated  from  the  oil  pressure  and 
the  model. 

The  creep  model  accounts  for  the  effects  of  the  air-jack-loading  rate.  Although  an  individual  test  could 
be  modeled  well,  the  standard  deviation  of  A,  B,  and  xG  for  the  eight  valid  tests  was  37%,  57%,  and 
64%  of  these  values,  respectively.  The  mean  values  chosen  are: 

A  =  2.43  B  =  -0.037  min."1  xQ  =  20  min.  (2) 

Using  these  mean  values,  the  standard  deviation  of  the  calculated  value  of  the  maximum  air  pressure 
during  the  eight  ramp  tests  is  64  kPa.  The  response  of  one  test  is  shown  in  Figure  1.  The  calculated 
time-dependent  response  of  the  stress  is  shown  together  with  the  applied  airjack  pressure  and  the 
reading  from  the  gauge.  Figure  1  illustrates  clearly  why  the  time-dependent  response  is  important 
because  the  maximum  applied  stress  is  three  times  the  value  read  by  the  gauge.  The  response  in  Figure 
1  is  typical  in  that  the  increasing  part  of  the  pressure  curve  is  matched  better  than  the  decreasing  part. 
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Figure  1.  Ramp  Test. 

A  theory  for  the  stresses  due  to  temperature  gradients  in  an  unconfined,  thick,  elastic  plate  of  uniform 
thickness  in  which  the  temperatures  and  the  properties  are  a  function  only  of  depth  is  given  by 
Timoshenko  and  Goodier  (1951): 

aEAT  1  h  2  12v  h  2 
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where  x  and  y  are  the  coordinates  in  the  plane  of  the  plate,  z  is  the  coordinate  down  from  the  mid-depth 
of  the  plate,  a  is  the  stress  at  z,  a  is  the  thermal  expansion  coefficient,  E  is  the  elastic  modulus,  AT  is 
the  temperature  change  from  a  uniform  temperature  as  a  function  of  z,  v  is  the  Poisson’s  ratio  for  the 
material,  and  h  is  the  plate  thickness.  To  determine  the  thermal  stress  on  the  stress  sensor,  the 
temperature  change  at  the  sensor  depth  was  used  to  evaluate  the  first  term  on  the  right  side  of  Equation 
(3).  For  sea- ice,  a  is  essentially  constant,  and  E  is  taken  to  be  zero  for  the  bottom  20%  of  the  depth, 
after  data  from  Cox  and  Weeks  (1988).  The  thermal  stress  in  sea  ice  is  thus 
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The  ice  temperatures  were  measured  with  thermistors  at  several  depths.  To  account  for  the  creep  out  of 
the  thermal  stress,  the  temperature  change  was  calculated  not  from  an  initial  uniform  temperature,  but 
from  a  reference  temperature  calculated  by  smoothing  the  measured  temperatures  with  a  low-pass 
filter.  The  filter  has  a  14-day  time  constant,  which  is  the  same  time  constant  as  the  slow  creep  found  in 
the  pressurization  decay.  Constant  values  of  1  GPa  for  the  effective  elastic  modulus,  0.3  for  Poisson’s 
ratio,  and  5.1  x  10"5  per  °C  for  the  thermal  expansion  coefficient  were  found  to  make  the  calculated 
thermal  stress,  which  is  the  same  in  all  (horizontal)  directions,  consistent  with  the  ice  stress  measured 


with  all  four  sensors  in  the  rosettes.  These  values  are  similar  to  those  used  by  Lewis  (1993)  in  his 
calculations  of  thermal  stress  in  sea  ice.  Lewis,  however,  uses  a  variable  elastic  modulus,  a  different 
definition  for  a  reference  temperature,  and  creep  strain  in  his  model.  The  result  of  this  calculation  is 
illustrated  in  Figure  2.  The  effect  of  ice  thickness  change  has  been  accounted  for.  The  ice  was  76  cm 
thick  on  Day  -35  and  115  cm  thick  on  Day  57. 


IMPACT  FOR  SCIENCE  AND/OR  SYSTEMS  APPLICATIONS 


The  anisotropic  sea-ice  mechanics  model  can  be  developed  for  application  to  a  range  of  problems,  such 
as  assisting  in  the  interpretation  of  SAR  data,  forecasting  ice  conditions,  using  current  SAR  data, 
forecasts  of  ice-generated  noise,  studies  of  climate  dynamics,  ice  loads  for  offshore  platform  design, 
and  environmental  hazard  analysis. 

TRANSITIONS 

At  present,  NWRA  is  using  the  sea-ice  stress  data  collection  and  interpreted  under  ONR  programs  for 
evaluating  the  new  ice  model. 

RELATED  PROJECTS 

The  project  is  being  done  together  with  a  project  with  R.  S.  Pritchard  of  IceCasting,  Inc.  All  of  the  data 
has  been  sent  to  him  for  his  use  in  model  development. 
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